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We study color fluctuations in the quark-gluon plasma produced at the early stage of nucleus- 
nucleus collision at RHIC or LHC. The fluctuating color current, which flows along the beam, can be 
very large due to the strong anisotropy of the parton momentum distribution. A specific fluctuation, 
which splits the parton system into the current filaments parallel to the beam direction, is argued 
to grow exponentially. The physical mechanism responsible for the phenomenon, which is known as 
a filamentation instability, is discussed. 

o\ ■ 

In the near future the nucleus- nucleus collisions will be studied experimentally at the accelerators of a new gen- 
eration: Relativistic Heavy-Ion Collider (RHIC) at Brookhaven and Large Hadron Collider (LHC) at CERN. The 
collision energy will be larger by one or even several orders of magnitude than that one of the currently operating 
machines. A copious production of partons, mainly gluons, due to hard and semihard processes is expected in the 
heavy- ion collisions at this new energy domain Q. Thus, one deals with the many-parton system at the early stage 
of the collision. The system is on average locally colorless but random fluctuations can break the neutrality. Since 
the system is initially far from equilibrium the color fluctuations can noticeably influence its evolution. 

(NJ , In our previous papers Q we have studied the plasma oscillations in the nonequilibrium quark-gluon system pro- 
duced in ultrarelativistic heavy-ion collisions. We have argued that due to the strong anisotropy of the parton 
momentum distribution there are unstable modes which exponentially grow in time. Here we show that the fluctua- 
tions which initiate these modes are large, much larger than in the equilibrium plasma. We also discuss the physical 
mechanism responsible for the growth of the fluctuation which splits the parton system into the color current filaments 

^\ • parallel to the beam direction. For completeness we recapitulate at the end of this letter the results from our papers 

"*^h 0, which show the existence of the unstable modes. 

The distribution functions of quarks Qij(t,x,p), antiquarks (t, x, p) , and gluons G a b(i,x,p) with i,j = 1,2,3 
and a,b = 1, 2, 8 are |^^] matrices in the color space; t, x and p denote the time, position, and the momentum. 
The color current expressed through these functions reads |||| 

«(<,*) = i» j if (»s(<3«(f,x,p)-Q iJ (<,x,p))+ir fc G l ,„((,x, P ) 

X I where , „ the QCD couphpg eons J, r- I the SU {3) s ™»p s »e„to r , /- the „pec.„e _ co-ant, and 
p M = (E p ,p) is the parton four-momentum with E p = |p| being the energy of the massless quark or gluon. 

We assume that the quark-gluon plasma is on average locally colorless, homogeneous, and stationary. Thus, the 
distribution functions averaged over ensemble are of the form 

(Qij(t,x,p)) = S ljn (p) , (Qij(t,x,p)) = <%n(p) , (Gij(t,x,p)) = %%(p) , 

which give the zero average color current. 

We find the fluctuations of the color current generalizing a well-known formula for the electric current For a 
system of nonintcracting quarks and gluons we have in the classical limit the following expression 

Mj(*,x) d ^ (J5(*i,xi)i?(* 2 ,x 2 )> = l?S*J^W /(p) (x- vt) , (1) 

where the effective parton distribution function f(p) equals n(p) + ^(p) + 6n s (p), (t, x) = (t 2 — £i,X2 — xi), and 
v = p/Ep is the parton velocity. Due to the average space-time homogeneity the correlation tensor depends only on 
the difference (t 2 — t\, x 2 — Xi). 
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The physical meaning of the formula ([!]) is transparent. The space-time points (ii,xi) and (t2,X2) are correlated 
in the system of noninteracting particles if the particles fly from (ii,xi) to (*2,X2). Therefore the delta <5( 3 )(x — vi) 
is present in the formula (|l]). The momentum integral of the distribution function simply represents the summation 
over particles. 

One finds the Fourier spectrum of the fluctuations from eq. (|l|) as 

M£>, k) = \g 2 6" h J ^ ?jL /(p ) 2nS(u kv) , (2) 

where u> is the frequency and k the wave- vector. 

We model the parton momentum distribution at the early stage of ultrarelativistic heavy-ion collision in two ways: 

f(p) = ^(X-y)e(Y + y)h(p ± )^— 1 (3) 
2Y p± chy 

and 

/(p) = 2^e(p-j>u)e(p+p||) Kp±) , (4) 

where y, pu, p±, y, and <fr denote the parton rapidity, the longitudinal and transverse momenta, and the azimuthal 
angle, respectively. The parton momentum distribution (Q) corresponds to the rapidity distribution which is flat in 
the interval (— Y, Y). The distribution (^) is flat for the longitudinal momentum — V < p» < V . We do not specify the 
transverse momentum distribution h(p±), which is assumed to be of the same shape for quarks and gluons, because 
it is sufficient for our considerations to demand that the distributions (||, |i|), are strongly elongated along the z— axis 
i.e. e Y > 1 and (p||) > (p±). 

Due to the symmetry /(p) = /(— p) of the distributions the correlation tensor M^ v is diagonal i.e. M' 1 " = 

for p ^ v. Since the average parton longitudinal momentum is much bigger than the transverse one, it obviously 
follows from eq. (^) that the largest fluctuating current appears along the z— axis. Therefore, we discuss the M zz 
component of the correlation tensor. M zz (ui, k) depends on the k— vector orientation and there are two generic cases: 
k = (fe x ,0, 0) and k = (0,0, k z ). The inspection of eq. (||) shows that the fluctuations with k = (fc^O, 0) are much 
larger than those with k = (0, 0, k z ). Thus, we compute M zz (ui, k x ). 

Substituting the distributions (|^,||) into (j^) one finds after azimuthal integration 

j \ rab g 2 f Y , f°° , s sh 2 y <d(k 2 x -u 2 ch 2 y) 

M%(uj,k x ) = S ab — \- / dy / d P± h(p±)px -rr 1 > 5 

32^ Y J_ Y J chy J k 2 2 ch 2 



f v f°° P±P\\ ®(k 2 p 2 t -u 2 El) 

One observes that the integrals from eqs. (HJf]) reach the maximal values for lu 2 <C k 2 . Therefore, we compue M z £ 
for uj = 0. Keeping in mind that e Y 3> 1 and (pii) S> (pi_) we get the following approximate expressions for the flat 



y— and pu— distributions: 



M*J(u; = 0,k x ) = l9 2 S ab j^- M- (8) 
8 (px) \kx\ 

where (p) is the effective parton density given as 

(p) = / /(p) = 'h 2 1 dp±p± - h( - p ^ = I ^ p)qq + 1 ' 

with (p)qq denoting the average density of quarks and antiquarks, and (p) g that of gluons. For the flat pu— case we 
have also used the approximate equality 
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I dp±_ h(p±) = - — r / dpj_p A 
Jo (P-l) Jo 



h{p a 



to get the expression (g) 

It is instructive to compare the results (f7|, ^J) with the analogous one for the equilibrium plasma which is 

One sees that the current fluctuations in the anisotropic plasma are amplified by the large factor which is e Y /Y or 

Let us now discuss how the fluctuation, which contributes to M z §(uj = 0, k x ), evolves in time. The form of the 
fluctuating current is 

ja(x) = ja e z cos(k x x) , (9) 

where e z is the unit vector in the z— direction. Thus, there are current filaments of the thickness 7r/|fc^| with the 
current flowing in the opposite directions in the neighboring filaments. 

In the limit of weak fields the chromodynamics can be approximately treated as an eight-fold electrodynamics. 
Consequently, the magnetic field generated by the current (^) is given as 

B Q (x) = j- e y sin(k x x) . 

fox 

The Lorentz force acting on the partons, which fly along the beam, equals 

F(x) = q a v x B a (x) = -q a v z ^- e x sm(k x x) , 

K-x 

where q a is the color charge. One observes, see Fig. 1, that the force distributes the partons in such a way that those, 
which positively contribute to the current in a given filament, are focused to the filament center while those, which 
negatively contribute, are moved to the neighboring one. Thus, the initial current is growing. 

The mechanism described here is well-known in the plasma physics [Q and it leads to the so-called filamentation 
or Weibel instability j6|. In the context of the quark-gluon plasma the phenomenon has been first discussed in the 
system of two interpenetrating parton beams ||. Such a system however seems to be completely unrealistic from 
the experimental point of view. Then it has been argued [^) that the filamentation instability can occur under 
the conditions which will be realized in heavy-ion collisions at RHIC and LHC. Let us briefly recapitulate these 
considerations which confirm the qualitative arguments presented above. 

The spectrum of plasma modes initiated by the color fluctuations is determined by the dispersion equation which 
for the anisotropic plasma is 

det|k 2 ^' - k l k j - cA y '(w,k)| = , i,j = x,y,z 
where the chromodielectric tensor e lJ is 

d 3 p v l df(p) 



kv + i0 + dp 1 



kv\ k l v^ 



LO 



For the specific color fluctuation with wave vector along the x— axis and the electric field parallel to the z— direction, 
the dispersion equation simplifies to 

H(uj) = k 2 x -uj 2 e zz (uj,k x ) = 0. (10) 

The so-called Penrose criterion states that the dispersion equation H(lu) = has unstable solutions if H(lu = 0) < 0. 
We have shown that under the reasonable assumptions concerning the form of the transverse momentum distribution 
h{p±), the criterion is indeed satisfied. Then, we have solved approximately the dispersion equation ( |To| ) and found the 
unstable mode with the pure imaginary frequency. The estimated characteristic time of the instability development, 
which appears as short as 0.3 — 0.4 fm/c, is significantly smaller than the life time of the plasma created in the nuclear 
collision. Therefore, the instability can indeed develop and play a significant role in the system dynamics. 



3 



One asks whether the color instabilities are detectable in ultrarelativistic heavy-ion collisions. The answer seems to 
be positive because the occurrence of the filamentation breaks the azimuthal symmetry of the system and hopefully 
will be visible in the final state. The azimuthal orientation of the wave vector will change from one collision to another 
while the instability growth will lead to the energy transport along this vector (the Poynting vector points in this 
direction). Consequently, one expects significant variation of the transverse energy as a function of the azimuthal 
angle. This expectation is qualitatively different than that based on the parton cascade simulations where the 
fluctuations are strongly damped due to the large number of uncorrelated partons. Due to the collective character 
of the filamentation instability the azimuthal symmetry will be presumably broken by a flow of large number of 
particles with relatively small transverse momenta. The jets produced in hard parton-parton interactions also break 
the azimuthal symmetry. In this case however the symmetry is broken due to a few particles with large transverse 
momentum. The problem obviously needs further studies but it seems that the event- by- event analysis of the nuclear 
collision give a chance to observe the color instabilities in the experiments planed at RHIC and LHC. 
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Figure captions 

FIG. 1. The mechanism of filamentation. The phenomenon is, for simplicity, considered in terms of the electrody- 
namics. The fluctuating current generates the magnetic field acting on the positively charged particles which in turn 
contribute to the current (see text). (g> and © denote the parallel and, respectively, antiparallel orientation of the 
magnetic field with respect to the y— axis. 
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